Bulk and edge correlations in the compressible half-filled quantum Hall state 
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We study bulk and edge correlations in the compressible half-filled state ji],^], using a modified 
version of the plasma analogy. The corresponding plasma has anomalously weak screening proper- 
ties, and as a consequence we find that the correlations along the edge do not decay algebraically 
as in the Laughlin (incompressible) case, while the bulk correlations decay in the same way. The 
results suggest that due to the strong coupling between charged modes on the edge and the neu- 
tral Fcrmions in the bulk, reflected by the weak screening in the plasma analogue, the (attractive) 
correlation hole is not well defined on the edge. Hence, the system there can be modeled as a free 
Fermi gas of electrons (with an appropriate boundary condition) . We finally comment on a possible 
scenario, in which the Laughlin-like dynamical edge correlations may nevertheless be realized. 
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I. INTRODUCTION AND PRINCIPAL RESULTS 

Laughlin's theory of the fractional quantum Hall effect 
(QHE) H was given in terms of wave functions of the 
ground state and quasihole excitation. Using a plasma 
analogy to calculate the static many-body correlators, 
which characterize these wave functions, he was able to 
advance a very successful physical picture of the electron 
system. The wave functions, describing the incompress- 
ible states, contain the Laughlin-Jastrow factor, which 
leads to a special, later introduced, Girvin-MacDonald 
(GM) correlations in the bulk Q , and Wen's correlations 
on the edge H|). The Laughlin-Jastrow factor is everp- 
resent in QHE states - it exists even in the compressible 
half-filled state jjj, for which an explicit wave-function 
has been proposed by Read and Rezayi (RR) Q). The 
question arises whether its manifestations, in terms of 
the above mentioned correlations, survive in more gen- 
eral quantum Hall states, and in particular in the com- 
pressible states. Why is this question important? The 
correlations that are embodied in the Laughlin-Jastrow 
factor lie at the heart of various quasiparticle pictures 
(composite fermions, composite bosons) of the QHE 
in the bulk. From the theoretical viewpoint, it is inter- 
esting to understand the status of Bose condensation, 
implicit in the Laughlin-Jastrow factor , in the com- 
pressible state. Related to this is the question to what 
extent Laughlin's quasihole construction in the compress- 
ible state (a zero of the wave function) can be considered 
as an elementary excitation of the system. 

Experimentally, these correlations are in principle ac- 
cessible by tunneling measurements. Indeed, recent 
edge-tunneling experiments by M. Grayson et al |lC| ] 
prompted the question whether the Luttinger liquid pic- 
ture |,(|, which is characterized by Wen's correlations, 
is valid for general quantum Hall systems, including 
the compressible states. A number of theoretical works 
ITLEliI have attempted to explain the puzzling results of 



Ref. |L0[ , in terms of charged excitations on the edge that 
are effectively decoupled from the bulk |13|. 

In this paper we concentrate on the compressible QHE 
system at filling factor one-half. We assume that the 
RR wave function well describes the ground state of the 
system, even when we consider a system with an edge. 
Namely, we assume the composite fermion (or more pre- 
cisely dipole) picture j7|,[| to apply everywhere. We 
rederive the GM and Wen's correlations in the Laugh- 
lin state considering the leading order contributions of a 
weak-coupling plasma approximation (see also Ref. Jill ) . 
Then we consider the same correlations (appropriately 
redefined) in the RR state. In calculating these we use 
the same approach - a systematic expansion of a plasma 
free energy - with necessary modifications to include the 
Fermi sea correlations This introduces a statisti- 

cal mechanics viewpoint of the problem, in terms of an 
anomalous, weakly-screening plasma. 

Applying the forementioned procedure (and view- 
point) on the RR state, we find that Wen's correlations of 
the edge do not decay algebraically (at large distances) as 
in the Laughlin state. This excludes the possibility of ex- 
istance of a subspace of charge density waves on the edge 
(of the type found in the Laughlin state) , that is decou- 
pled from the rest of the excitations - i.e., the neutral 
bulk excitations. The form of the obtained equal-time 
electron Green's function on the edge suggests that, in 
the first approximation, the physical picture of the RR 
edge is that of a Fermi gas of electrons. The bulk GM 
correlations, on the other hand, decay algebraically, in 
an almost identical way as in the Laughlin state. 

Below we detail the derivation of the correlators, in the 
bulk (Sec. [n]) and on the edge (Sec. III). A discussion of 
theoretical and experimental implications of the results 
is given in Sec. [FvL 
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II. CORRELATIONS OF THE BULK 

In this section, we employ the plasma analogy to de- 
rive the appropriately generalized GM correlator in the 
compressible RR state. To introduce the method, we first 
use it to derive the known result for the Laughlin state 
(Eq. (|) below). 

A. Correlations of the bulk in the Laughlin state 

In the Laughlin state, corresponding to filling factors 

— with m odd, the GM correlator M is defined as the 
density matrix, 

p(z,z') = 

N j d 2 z 2 ■■■ J d 2 z N ^ b (z, z 2 ,..., z N ) x Vbjz', z 2l ■ ■ ■ , z N ) 
Jd 2 Zl ---fd 2 z N \^\ 2 ' 

(1) 

for the bosonic many-body function, 

*5 = II 1^ - ^r i exp{-i^ |^| 2 }, (2 ) 

i<j 

obtained from the Laughlin wave function by omitting 
the phases of the relative distances between any two elec- 
trons, (zi — Zj). As shown in Q|, the asymptotic form of 
p(z, z') is 

p(z,z')~\ z -z'\-"- (3) 

This correlator expresses a Bose condensation, with alge- 
braic off-diagonal long range order, of composite bosons 

- defined as electrons with m flux quanta attached. We 
now derive the above form using the weak-couplig plasma 
analogy. 

We first rewrite the integrand as 0] 

® b (z,...,z N ) x ^ b (z',...,z N ) = 
i 

exp{2m In \zi — Zj\} 

Kj 

i i 
x exp{+m In \ z — Zi\ + m^~^ ln\z — z,-|} 

i i 

^e M - l -j2\z^-\\zr-\\z?} 

(4) 

and similarly the numerator. (The prime means that i=l 
is excluded from the summations.) Using the Laughlin 
plasma analogy we can write p(z, z') as 

P (v') = k^r f S4», (5) 

Zi(Z,Z) 



where Z(z, z') is a partition function of a classical 2D 
plasma with inverse temperature /3 = — , each particle 
with charge m, and two impurities with charge ^ each, 
at the locations z and z' . (Z(z, z') is a partition function 
with one impurity of charge m at an arbitrary location, 
because the value of the partition function does not de- 
pend on z.) n is the average density of particles (equal 
to tt— in the usual units). To calculate the ratio of the 
two partition functions we may expand the exponentials 
in the parameter m, which we will assume to be small. 
The expansion will generate terms that can be described 
by diagrams and corresponding rules. 

As usual in this kind of expansion in the statistical 
mechanics analogue, the expansion of the denominator 
involves only connected diagrams. Each diagram consists 
of parts, hereon called disconnected parts, which connect 
two impurities at z and z' but are otherwise disconnected 
among themselves. Then, the rules that correspond to 
each diagram in the expansion are as follows: 

(1) Associate with each interaction line a two- 
momentum satisfying momentum conservation at each 
internal vertex, 

(2) Associate with each interaction line between par- 
ticles — 27r %™ , with each interaction line between a par- 

|fe| 

ticle and an impurity — 27r ^ 2 X 2 , with each interaction 

2 *K Q ( — 

line between impurites -4j , and with each internal 

\k\ 

vertex n, 

(3) For each incoming (from z) (which is also outgoing 
to z') momentum for each disconnected part, integrate as 
J h^ts exp{ifc(r — r*)}, but for each internal momentum 

(4) Multiply with a symmetry factor (if any). The sym- 
metry factor is an inverse of the number of ways that we 
can interchange a given number of identical parts of a 
given diagram, and recover the same graph. 

The diagrams that represent the ineraction with the 
background are mutually canceled (as we checked for the 
first diagrams in the expansion) and we will not consider 
them. In our problem the density n is fixed and depends 
on the small parameter m. In order to get the correct 
order of the diagram (i.e. the power of m) in the expan- 
sion, we must take this into account. The lowest order 
diagram has value one. The next in order are diagrams 
of the form shown in Fig. 1 and are of order m. We can 
easily sum them and the result is 

V,,,{\?-i'\) = 

|/c| 2 

The sum represents an effective screened interaction be- 
tween two impurities. The infinite summation of cer- 
tain type of diagrams that diverge even more singulary 
as we increase the number of interaction lines, is a well 
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known ansatz in the many-body theory of the Coulomb- 
interacting electron gas in three dimensions. This cap- 
tures well the phenomenon of screening that is charac- 
teristic of long-range forces. In our case the infinite sum- 
mation is even further enforced, given the fact that the 
diverging diagrams are of the same order in m. 



FIG. 1. The diagrams leading to screening of the interac- 
tion in the bulk. 

We now rewrite the ratio of partition functions in Eq. 
(i as 



Zi{z, z) 



(7) 



where Af(z,z') represents the difference in the free en- 
ergy between the two configurations of the impurities. 
The above exponential form can be obtained by sum- 
ming the set of diagrams, whose disconnected parts are 
of the form shown in Fig. 1. Thus we find 



Z(z,z') 
Z(z, z) 



= exp{V e// (r-r')}. 



(8) 



As 



oo the ratio approaches unity, because V e 

Hence, 



is an effective screened interaction |16 



Eq. 



reduces to the well-known expression for the GM corre- 
lator, Eq. (g). 

This result is derived, and found to have the same 
form for larger, physical m's [ij. Therefore, it is pos- 
sible to analytically continue the correlator obtained in 
the weak-coupling approach to larger m's. Applying the 
same weak-coupling infinite summation, it can be shown 
that the continuation is valid also in the calculation of 
the static structure factor in the small-momentum limit 
(when corrections to the infinite summation are added, 
this includes also the term proportional to the fourth 
power of the momentum [fhH ) . 

It is interesting to check what does the weak-coupling 
approach yield for the distribution of the charge in the 
tail of the Laughlin quasihole excitation || . The quantity 
that describes this is 101 



gv2\\zi - W\ 



N f f N 

— J d 2 z 2 ---J d 2 z N \ J[(zi - tu)* 



(9) 



where "J and Z is the Laughlin wave function and its 
norm, respectively. In order to capture the physics of 



screening we sum the same, most imortant diagrams as 
before, and approximate Eq. (g) by [|l4| 



gu(\zi - w\) = n 
d 2 k 



exp{ifc(ri — w)}- 



2nP(mxl) 
\k\ 2 



1 + n 



(10) 



As \z\ — w\ — > oo the function 312 — w\) should tend 
to the unpertubed density n, and it behaves as 



gi2(\zi - to I J = n 
1 



-m Const- 



V\zi - w\ 



exp{ } 



(11) 



where {\/r 2 D ) = 2ir(3m 2 n = 2, m being the Debye length. 



B. Correlations of the bulk in the compressible 
half-filled state 



The theory and physical picture of the filling frac- 
tions — where m is even, evolved from some Fermi con- 

m 5 

densation of charged (Chern-Simons) composite fermions 
(electrons with even number of flux quanta attached) to 
a well-defined Fermi condensation of dipole quasiparti- 
cles p|. This emphasized the advantage of Read's pic- 
ture ]17[, which, from the begining, takes into account 
the binding of electrons to (so-called) correlation hole. 
(Equivalently, the statement is that the zeros of the many 
body functions are found at or near the electrons). At 
even denominators the overall neutral composite object 
is a dipole (with Fermi statistics). 

The ground state wave function that corresponds to 
this picture is the RR wave function || 



= 7 ? LLi{det[exp(i£; l i? i )]1' L }, 

h3 



(12) 



with a Slater determinant of free waves that fill a Fermi 
sea, which when projected to the lowest Landau level 
(LLL) (Vlll stands for the projector), acts on - 
the Laughlin wave function. In Eq. (|l^) we wrote 
the determinant in terms of plane waves, which con- 
stitute a convenient basis for a system of free parti- 
cles (in a rectangular geometry). The Laughlin wave 
function, on the other hand, is very often expressed 
in the rotationaly-symmetric gauge (corresponding to 
a rotationally-symmetric geometry) amenable to the 
Laughlin plasma analogy. In order to facilitate our 
computations we will keep these two distinct geometry 
choices in the RR wave function. We justify this by the 
fact that, first, we will be interested in the (longwave- 
length) properties of the system in the thermodynamic 
limit (when the boundary conditions should not matter), 
and second, each component of ^rr will enter our calcu- 
lations in the form of translationally invariant, geometry 
independent elements. 
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To illustrate [[2|J9[] the dipole physics contained in Eq. 
(|l2l), we note that the LLL projection translates factors 
of the form exp{i(fc"z)/2}, where k = k x +ik v and z = x + 
iy, into the shift operator exp{ikJ^}, which acts on the 
original (before projection) holomorphic (z - dependent) 
part of the wave function. This effectively means that 
each electron becomes displaced from the position of its 
correlation hole by (— ik) (where k takes values from the 
Fermi sea), and therefore dipole moments are induced. 

In the calculation of correlation functions the effects of 
the LLL projection can be taken into account by using 
the following identity 



/ 



d 2 zexp{— -|z| 2 } exp{iqr} 



d z* 
x {exp{-ifcj — } exp{-i ki — }F 1 (z*)} 
oz* 2 

d z 
x{exp{ik 2 — }exp{ik 2 -}F 2 (z)} 

r k*q, , k 2 q* k 2 k* 
= exp{ — }exp{+-^-}exp{ — } 

d 2 z exp{— — |z| 2 } exp{iqr} exp{— ikir\ 
exp{ik 2 r}F 1 (z*)F 2 (z) . 



(13) 



If we search for the long distance behavior of the corre- 
lation functions, usually the calculations give the same 
result as obtained from the unprojected version of the 
RR function. 

This is the case with the appropriately generalized GM 
correlations to the compressible case. The many-body 
wave function employed in the calculation of the density 
matrix Eq. ([!]) is 

$>(z, z 2 , . . . , z N ) = 

N 



X 



JJ \zi - Zj + ik a{l) - ik^j) | 2 
Y[\z- z t -ik a{l) \ 2 exp{-jY,\ z *\ 2 } ■ ( 14 ) 



N 

x n 1 - ~ ■' ' 

i=2 



We now introduced a particle, with coordinate z, with- 
out the (projected) plane wave that enters the Fermi sea 
part, therefore without the fermionic statistics that char- 
acterizes the rest of the (N — 1) particles. The rest of its 
correlations with other particles is as of any other parti- 
cle. Similarly to the Laughlin case, the phase part of the 
Jastrow-Laughlin factor with coordinate shifts is omit- 
ted. This, in the Chern-Simons picture, corresponds to 
attaching of two flux quanta (at distance (ik)) to each 
electron. 

Nevertheless, as can be shown, for the type of calcu- 
lations that we do, the projection to the LLL does not 



affect the final result and, for the sake of simplicity, we 
will explain the method on the unprojected version for 
which 

%(z, z 2 , ■ ■ ■ , zn) = 

N - 1 

JJ \z - Zi\ 2 det{exp(ikifj)} J[ \z t - Zj \ 2 exp{-- ^ \z t \ 2 }. 

i 

(15) 



i<3 



We next assume that the dominant correlations lie in 
the (Jastrow-Laughlin) differences, and for the moment 
neglect the Slater determinant. The complete plasma 
analogy is again possible and, as explained above, the in- 
finite summation of the diagrams of type Fig. 1 (for small 
to) is relevant. In the presence of the determinant the 
first necessary correction to this picture is the introduc- 
tion of a new vertex, that captures also possible Fermi 
(exchange) correlations between two points in the coor- 
dinate space |l4j| . In the momentum space, the vertex 
then corresponds to the static structure factor of the free 
Fermi gas, 



so(q) 



d 2 r exp{iq ■ r}(g(\r\) 



1) 



(16) 



where q 7^ 0, and the radial distribution function is 
9(r) = 

1 f d 2 h f d 2 k 2 t t 

n JktEF.S. ( Z7r ) Jk 2 £F.S. ( Z7T ) 

(17) 

(F.S. stands for the Fermi sphere). Symbolically, the new 
vertex is depicted in Fig. 2, as a sum of a direct and an 
exchange part, in which full lines represent Fermi particle 
lines. 



= o + (X^o 

FIG. 2. The vertex so(q) in the bulk of a RR state. 
From Eq. (|l7|), and the definition Eq. (|l6|), 



s (q) - n 



d 2 k 



(18) 



Here R represents the area of overlap between two Fermi 
spheres as shown in Fig. 3, where the center of one of 
the two spheres is displaced by q from the center of the 
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other one. The value of So(q) is then given exactly by 
the shaded area in Fig. 3. The area is easily calculated 
for \q\ small and the result is 



3k f \q\ 



(19) 



With the necessary introduction of the new vertex so(q), 
the interaction becomes less effectively screened than in 
the usual (Laughlin) case. It becomes 



VeffM) 



1 



27T/3TO" 
W 



so(q) 



(20) 



where the denominator can be interpreted as an anoma- 
lous dielectric constant of the corresponding modified 
plasma. In the coordinate space, at large distances, 
V e ff ~ 1/r, i.e. it is still long ranged and only partially 
screened. 




no available results for larger m to compare with. It was 
found 1 14 1 that the weakly-screening plasma approach 
gives the right (valid also for large m ]1S|] ) leading (small- 
momentum) behavior for the static structure factor of the 
compressible state, and generates expected (odd) pow- 
ers of momentum in the expansion. If we try to go be- 
yond the approach, and look for small- m (expected) cor- 
rections, it seems that they can not be generated [ fl4) . 
This is probably due to the nonanalyticities present in 
the compressible case (which were absent in the Laugh- 
lin case) that do not allow a perturbative treatment. 
Therefore we believe that the approach (essentially non- 
perturbative) can generate the correct (large m) lead- 
ing behavior for the correlations that we study. They 
are between points which are directly connected only to 
the charge (Jastrow-Laughlin) part of the wave function, 
and that immediately suggests an approach that captures 
screening for their calculation. 

The weak-screening property of this modified plasma 
can be very well seen by considering a zero of the electron 
coordinates at a point w 0], which corresponds to the 
Laughlin quasihole in the incompressible case, 



(21) 



To simplify the calculation of the distribution of charge 
in the tail of this excited state, we will assume the un- 
projected version of the RR state in Eq. (§) (The use 
of the projected state involves some complications which 
are not essential and do not influence the final result). 
Now the appropriate infinite sum of the modified plasma 
can be expressed, in terms of diagrams depicted in Fig. 
4, with the shaded circle representing the new vertex, i.e., 
in this case, 



ffia(|«i -w\) = n 
d 2 k 



exp{ifc(ri — ?«)}■ 



1 



*o(fc)^|r 



(22) 



FIG. 3. The overlap of two shifted Fermi spheres. 

Nevertheless, if (keeping this change in mind) we apply 
the same summations and arguments as in the Laugh- 
lin case, we come up with the same algebraic decay of 
the GM correlations as in that case p4| . This decay 
is slightly modified by the exponential (Eq. (||)) of the 
partially-screened interaction (effectively a constant as in 
the Laughlin case at large distances). 

The question that arises immediately is whether the 
analytical continuation to larger (physical) m's is pos- 
sible, and, moreover, whether the screening plasma ap- 
proach is reliable in giving the leading behavior of the 
correlator. Because of the absence of a complete anal- 
ogy with some physical, well-studied plasma, there are 



FIG. 4. The diagrams contributing to the calculation of 
512 (|^i — w\) in the RR state. 

The most important contributions to g±2 in the limit 
\zi — w\ — > oo come from nonanalyticities present in the 
integrand. They stem from the nonanalytic behavior of 
so at k = and k = 2kf. Assuming that the small- 
momentum result Eq. (|l9|) for sq is valid for any k 
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(analogously to the Thomas- Fermi approximation for the 
electron gas in three dimensions) we get the contribution 
from the k — region, 



and 



{512(^1 - w\) ~ n}\r.F. oc 



m \z\ — w\' c 



(23) 



The contribution from the k = 2kf region can be calcu- 
lated to be 

{512(^1 - w\) - n}\ F .o. oc -j--, — pr sin(2A:/r) . 

kf \zi — w[ } 

(24) 

We may conclude from the expressions Eq.(p3|) and Eq. 
(|J), which summed up give the change in the distribu- 
tion of the charge from the uniform ground state contri- 
bution n (in the \z\ — w\ — > 00 limit), that the density 
far from the point w tends to n very slowly in compar- 
ison with the Laughlin case. The charge of this excited 
state (which may be argued to be ^ as for the Laughlin 
quasihole | fl7|| ) is spread over a much larger region than 
the one in the Laughlin case, due to the poor-screening 
properties of the modified plasma. 



III. CORRELATIONS OF THE EDGE 
A. Edge correlations in the Laughlin case 

In ||, Wen showed how calculation of the equal-time 
correlator along the edge in the Laughlin case can be 
reduced to the problem of finding the electrostatic en- 
ergy of placing an impurity outside the Laughlin plasma. 
For the sake of completeness and easy reference for our 
calculation we will, in brief, repeat his arguments. We 
then demonstrate, that the result is recovered in a weak 
coupling expansion. 

Review of Wen's procedure. We consider a disc 
of the Laughlin plasma, with a fixed radius R, at fixed 
filling factor — . As we increase the number of particles 
N, the density will increase (with appropriate change in 
the magnetic field B to keep — constant) and the de- 
scription that neglects detailes of the order of a mag- 
netic lenght would be more and more accurate, and the 
Laughlin plasma will behave as a metal (with its screen- 
ing properties). 

To calculate the edge correlator we envision placing an 
impurity of charge m outside the disc of such a plasma, 
at a distance z where \z\ ^> R (so that the details of 
the edge do not matter), and consider the ratio Zi/Z, in 
which 

Zj(z,~z) — / "|~[ d 2 Zj exp{^~~^ 2m In \zj — Zj|}x 

- 1 

x exp{+^{--|z fc | 2 + 2mln \z - z k \}} (25) 



Z= f nd 2 z l exp{^2mln|z 4 -z J | +^{-i|z fc | 2 }}. 

(26) 



k=l 



From the first-quantization (quantum-mechanical) point 
of view the ratio is the one-particle (electron) density at 
point z. On the other hand, from the point of view of 
the plasma analogue, In -± is the electrostatic energy re- 
quired to transfer the impurity from infinity to the point 
z. This energy can be expressed as 



In- 



Z 



mN2 In \z\ — mln(l 



R 2 



)+o(l/N). (27) 



The first contribution is the electrostatic energy between 
the total charge N and the impurity, where in the first ap- 
proximation the plasma droplet is assumed undeformed 
by the presence of the impurity. The second contribu- 
tion describes the most important part of the deforma- 
tion that occurs: the image charges of the impurity [ p^[ . 
The rest of the contributions are expected to be of or- 
der jj or less (Due to the form of the first contributions, 
analiticity in TV is expected.) 

To find out the electron correlator between points z\ 
and Z2 (on the edge), Wen first noticed that the expres- 
sion on the right-hand side of Eq. (|27]) is holomorfic in z 
and anti- holomorfic in z (outside the system), and there- 
fore can be analytically continued, i.e. 



In — 1 *' — — w mNhi{z\Z2] — m\n{l 



Rr , 
Z1Z2' 



(28) 



z\ and Z2 can be considered to be even on the edge if 
the final result of the analytical continuation exists, i.e. 
if it is finite. This excludes the points Z\ — Z2 on the 
edge (|^x| — I #2 1 — K)i where the above expression is 
logarithmically singular. Then, if z\ — i?exp{i-^} and 
Z2 = R, the electron correlator is (in the disc geometry) 



< L\&{z x )y{z 2 )\L >= 

Zl{zi,Z 2 ) , 1, |2 , r 1, ,2-1 

7? exPl-jN }exp{--|z 2 | } 



(29) 



In the limit ^ <C 1, where circular and rectangular ge- 



ometries are indistinguishable, this becomes 



< L|#t(zi)¥(*jj)|L >r 



1 



(30) 



k=l 



which coincides with the correlations on the edge ob- 
tained in the (more familiar) bosonization approach. 

Derivation of the electrostatic energy of the 
edge impurity in the Laughlin case using the 
weak-coupling plasma expansion. According to 
Wen's idea, in order to find the equal-time electron corre- 
lator, it is sufficient to compute the electrostatic energy of 
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an impurity of charge m at a point z outside the Laugh- 
lin plasma. We now describe the diagramatic solution 
of this Statistical Mechanics problem. To simplify the 
calculation, we consider a plasma which extends over the 
half-plane x < instead of a disc (in the thermodynamic 
limit, the choice of geometry is immaterial); the impu- 
rity coordinate is z — £ (along the positive x-axis). The 
derivation of this electrostatic energy, using the weak- 
coupling plasma expansion, parallels that of the density 
in the bulk; i.e. calculating the electrostatic energy of 
a particle interacting with a negative background - the 
rest of the particles |20|. In the present case the system 
is not infinite in the x-direction, and that introduces a 
new type of vertex in the diagrammatic expansion. The 
vertex connecting two interaction lines of momenta qf 
in the x-direction, which in the infinite case is 



nS(qi - qf) 



(31) 



(where n is the density), is replaced in the half-plane case 
by 



2vr -z{qi-qf) 



(32) 



i.e. proportional to the Fourier transform of theta func- 
tion, 



exp{i(qi — qf)x}dx 



0(—x) exp{i(gi — qf)x}dx 



w in the half-plane are integrated over, corresponds to 
the first (direct term) in Wen's expansion, Eq. (|27|). It 
is also proportional to the size of the system, and strictly 
speaking diverges in the case of half-plane system. (This 
divergence does not matter, and can be handled by con- 
sidering a rectangular system with sizes L x and L y much 
longer than the distance (£) of the impurity from the y - 
axis.) 

The diagrams of the type depicted in Fig. 5(b) are 
not included (by using a screened instead of the bare in- 
teraction in Fig. 5(a)), although they are of order m 
as well. As we remarked earlier in this section, the dia- 
grams that should be taken into account are of the same 
form as the ones that we select to play the role of pos- 
itive background (i.e. those that cure divrgences in the 
expansion with the two-particle interaction) in the infi- 
nite system case. In that case, the diagram of the form 
in Fig. 5(a) cancels all divergences when the interaction 
line does not connect to any other interaction line. When 
the proper selection is done, and all diagrams that cure 
divergences are present, the complete partition function 
is well-defined and a constant. Similary with impuri- 
ties and in the semi-infinite case, if all due interactions 
(additional diagrams) are included in the partition func- 
tion (including the interaction of impurities with positive 
background) it becomes a constant (due to the screening 
property of plasma). The partition function Zj in Wen's 
derivation is not complete, and therefore the part on the 
right-hand side of Eq. ( p7j ) is not a constant, and can 
be associated with the interaction of the impurity with 
"negative background" . 



-i(Qi ~ If) 



+ 7r% 4 -q f ) 



(33) 



(a) 




a,) 



% (w) 5 (w) 



FIG. 5. (a) The leading contribution to the electrostatic 
interaction of a charged impurity with the plasma, (b) Dia- 
grams that are not included in Eq. (p7j). 

The diagrams that are leading in the small m expan- 
sion, and are of order m, are given in Fig. 5(a) and 
Fig. 6. The diagram in Fig. 5(a), in which the points 



FIG. 6. Diagrammatic representation of the image charge 
energy term. The vertex in the half-plane case (corresponding 
to Eq. ([32])) is denoted by a dotted circle. 

The diagrams in Fig. 6 are all relevant and deserve 
a special attention. Their value (at least in the long- 
distance limit), can be calculated by solving an integral 
equation for an effective vertex, V{qi,qf), 



V{qi,qf) = 

±{[ l - 

2vr -i(qi - q f ) 



dk[— 



1 



i(q* - k) 



+ TrS(q l - q f )] + 
+ ir5(qi - k)] 



-Anm 



(q 2 + fe 2 ) 



V(k,q f )}. (34) 
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This equation can be schematically introduced as in Fig. 
7, where we denoted only momenta along the x direction. 
The momentum q along the y direction is the same on 
every line as in the infinite-plane case. Then the con- 
tribution of all diagrams in Fig. 6, summarized by the 
diagram on the left hand side of Fig. 7, can be expressed 
as 



dq f dq 



(27T) 



(2tt) 



dqf exp{- 
-47rm 



(q 2 + qf) 



(q 2 + qj) 



<?/)£} x 
(35) 



The solution to the equation, given in the long-distance 
approximation, can be found in Appendix A. It repro- 
duces the electrostatic energy of the imputity and its 
image charge in the half-plane case, corresponding to the 
leading contribution to the second term in Eq. ( p7| ) when 
the disc is considered to be large (R ^> (\z\ — R)). 




FIG. 7. The infinite sum of diagrams included in Fig. 
6, represented as an integral equation for the effective vertex 
V(qi,qf) (dotted double-circle). 



B. Edge correlations in the compressible half- filled 
state 

Now we switch to the calculation of the edge correla- 
tions in the RR case using the diagrammatic method. We 
first consider the unprojected RR state, and as a basis of 
free waves that enter the Slater determinant, we choose 



exp(ik y y) cos(k x x) 



exp(ik y y) sin(k x x) 
V2n 



(36) 



where k x and k y take values from a Fermi box (not 
sphere) in the k - space. As in the Laughlin case, 
we assume that the radius of the Laughlin disc is very 
large in comparison with the distance (along the edge) 
over which we measure correlations. So, effectively, we 
again consider the half-plane problem for which, on the 



other hand, the basis choices (Eq.(|3q)) are also appro- 
priate; there is no discrepancy between geometries of 
the Laughlin-Jastrow and free-wave part in the ground 
state, as in the full-plane case. In Eq. ( ^6[) the coor- 
dinate x is measured from the edge of the half-plane, 
i.e. a tangent to the large disc. If, somehow, the charge 
and neutral (fcrmionic) part decouple on the edge, the 
choices (Eq. (p6|)) are quite natural, because they sat- 
isfy the requirement that the (neutral) current normal to 
the boundary is zero i.e., that the fermionic number is 
conserved. 

First, we consider the correlations of the object intro- 
duced in Eq. ( pi]) to which, due to the correspondence of 
its construction to the one of the Laughlin quasihole, we 
will refer as a quasihole. This, of course, does not entail 
that the quasihole is a well-defined object - eigenstate 
of the Hamiltonian, as in the case of the Laughlin quasi- 
hole. It might be such (on the edge) if we find that its 
correlations arc of the same type as in the Laughlin case 
(Eq. J3C|)), and, therefore the charge degrees of freedom 
(on the edge) in the RR state can be described in the 
Luttinger liquid framework (or, microscopically, by the 
possible states of quasiholes). Again, as in section III A, 
to mimic the charge part of the electron, in this case, 
we consider the correlations of the object constructed by 
putting m (m = 2) quasiholes at the same place. Then, 
to find out if there is a departure from the Laughlin case, 
we consider the density of this object outside the half- 
plane system described by the RR state. 

In the language of the modified plasma, we are placing 
a charged impurity (not directly connected to the plane- 
waves part) outside the system, and checking whether 
the image charge physics still holds. Due to the poor 
screening in this modified plasma, the charge induced 
by the external impurity does not accumulate near the 
edge (within a microscopic screening length) as in the 
ideal plasma. Rather, the induced charge is expected to 
slowly decay towards the interior of the system. To get 
a handle on the form of the electrostatic energy associ- 
ated with this effect, we can employ the Thomas- Fermi 
approximation to compute the induced charge, given the 
dielectric properties of the modified plasma derived in 
the previous section (Eq. ( ^0| ) and the prodeeding dis- 
cussion). The calculation is summarized in Appendix B. 
We find that unlike the Laughlin (ideal) plasma case, the 
leading behavior of the "image charge" electrostatic en- 
ergy is a constant, rather than a logarithmically singular 
term. We now derive this result systematically in the 
diagramatic expansion framework. 

We first must find an effective vertex which corre- 
sponds to Eq. ( j3^ ) in the Laughlin case, and represented 
by the diagrams in Fig. 6. That will be done at the same 
level of approximations as in the case of the bulk correla- 
tions. Explicitly, to find the value of the effective vertex, 
we consider the two contributions, direct and exchange, 
depicted symbolically in Fig. 8, in the simplest diagram 
with only two interaction lines. (The use of the dotted 
lines is to emphasize that we are now in the half-plane 
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FIG. 8. The modified half-plane vertex in the RR state; 
the second diagram represents the exchange contribution. 

The direct contribution (unapproximated) for the first 
choice of basis for the Fermi sea in Eq. (B6f) is 



= - I d A Zl \n z 



z x \{2m) n — 



1 E 

keF.S. 



cos 2 (k x xi) = 



= l(^) 2 J^J^x 
exp{-i(q x - p x )z}(2n)8(q y - p y ) x 



/ dxi exp{i(q x - Px)xi}j r J2k £ F.s. cos 2 (k x x 



(37) 



where the summation in k runs over the Fermi box. This 
summation can be rewritten as 



— cos 2 (k x x t ] 



keF.S. 



1 E 

keF.S. 



cos(2k x x\) 



(38) 



The second part can be neglected, because it leads to an 
effective smearing of both the 8 - function S(p x + q x ), 
and the pole ~ l/(p x + Qx) that we would get if it were 
a constant. The first part is the most important and sin- 
gular in the infrared limit, which dominates the infinite 
summation above. Therefore, the direct contribution to 
the effective vertex is 



1 n 



1 



2 27r -i(q x -p x ) 



+ ir8{q x -p x )} 



(39) 



i.e., half of the vertex in the Laughlin case. 

Applying similar arguments, that is keeping the most 
important terms that contribute to the value of the di- 
agram in the long-distance limit, we find that the ex- 
change contribution to the effective vertex is 



Qx)] 



(40) 



We get the same contributions, Eq. (39) and Eq. (40), 
for the second choice of basis in Eq. (ft6|). Therefore, the 
effective vertex in the RR case that parallels Eq. (^2|) in 
the Laughlin case, is 



—{ — 

2 27T l -z(q : 



1 



Px) 



Qx)} 



(41) 



It is not a simple multiple of the Laughlin vertex; be- 
cause of the exchange contribution, the solution of a new 
integral equation, corresponding to Eq. ( |34|) with the 
new "bare" vertex, will not yield the leading, logarith- 
mic behavior characteristic of the Laughlin case, which 
can be translated into the algebraic decay of the quasi- 
hole correlator. Namely, if a ^ 1 (for the definition of a 
see Appendix A), the solution of the integral equation in 
the long-distance limit is 



V( qi ,q f ) 



HQi - Qf) + 



^TSW^wil 2 + Qt)(Q 2 +Qf)- 



(42) 



The con- 



In our case b = — 47rm, a = i, and c = 4tt-. 

1 2 1 2 2n 

tribution from the 8- function to the electrostatic energy 
is 



mm — , 

Qc 



(43) 



where q c is the infrared cut-off. There is no obvious way 
to cancel the g c -dependence; i.e., in the case of the mod- 
ified plasma, we must keep the size of the system finite, 
and the distance of the impurity should be considered 
smaller than the size of the system, in the calculations. 
(Note that this appear to indicate, that the expansion in 
N is not analytic as in the Laughlin case.) The contri- 
bution of the second term can be written as 

+ m Const x (A - q c )q c x (g c f) + /(A, q c ) + o((q c £) 2 ) 

(44) 

where Const > 0, and / is an algebraic function of A and 
q c . In principle, other contributions, of order higher in 
m, from a more detailed solution of the integral equation 
can be calculated. We expect that their dependence on £ 
will be of the form (£,q c ) n or (j) n , where n takes on pos- 
itive integral values, (£ < ^-), and will not change the 
leading behavior Eq. (^3|) (in which we are interested the 
most). In the scope of our approach, which takes m small 
(and assumes the possibility of an analytical continuation 
to higher to), it is hard to estimate the true coefficients 
in front of the powers of £, due to the requirement to 
know them to all orders in to. Also, there might be rele- 
vant contributions from other diagrams (in the small to 
expansion) which we did not consider. But as we as- 
sume that the plasma correlations (although modified) 
are dominant for the calculation of the quasihole corre- 
lator, we do not expect that there will be any change in 
the leading behavior described by Eq. (|43|). 

The above calculations imply that the overlap between 
two quasihole excitations on the edge does not depend on 
the distance between them; it is a constant, but decreases 
with the size of the system. This might be understood, 
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taking into account that the quasiholes in the RR state 
are not well-defined, well-localized objects in the bulk 
(see the end of section II B), and certainly not on the 
edge where the screening seems to be even weaker than 
in the bulk. 

Once we take this point of view that, in fact, the states 
described by the Laughlin quasihole construction on the 
edge are extended, a special care must be taken concern- 
ing their normalization. In general, the normalization is 
expected to depend on the size of the system (as in the 
case of the free waves (in the non-interacting system)). 
Therefore, the first contribution to the plasma electro- 
static energy (Eq. ([13])) (that through the infrared cut- 
off depends on the size of the system) might be a conse- 
quence of an incomplete normalization of the quantum- 
mechanical correlator at the beginning of our calculation. 
If this term is included (in the normalization) from the 
beginning, the value of correlator at l arg e distances (in 
our approximation) approaches unity |21| . 

To find out the electron correlator, we must take into 
account the correlations that come from the neutral 
(plane-wave) part of the RR function, alone. These are 
not included in the preceding (modified plasma) calcula- 
tions, which gave the correlator of the quasihole, the ob- 
ject that (in our approximation) carries the charge part 
of the electron. The neutral contribution is expected to 
be of the form 



where g is a coefficient that depends on the boundary 
conditions. When combined with the charge correlator, 
it produces the usual (physical) decay of the electron cor- 
relator with the distance. Except for the dependence on 
the size of the system, the electron correlations on the 
edge are as if the system was a free (two-dimensional) 
Fermi gas of electrons. It can be shown that the same 
long-distance behavior of the correlations follows from 
the LLL projected wave function (Eq. (fl2|)). 

IV. DISCUSSION AND CONCLUSIONS 

If we assume that, indeed, the whole description of the 
edge of the compressible state is equivalent to that of 
a free Fermi gas, we can try to predict the occupation 
numbers (probability density) of electron near the edge. 
Then the second choice for the boundary condition in Eq. 
( |36| ) is more appropriate because the probability density 
(p{x)) should vanish at some point near the edge (x ~ 0). 
The resulting probability-density distribution 

, s r, sin(2fcpa;) 1 , . . 

p(?)~[kF ( ;C<0 )' ( 46 ) 

is very similar to the smooth function that one can get 
extrapolating the data that describe the occupation num- 
bers for electron near the edge in (finite-system) exact- 
diagonalization studies p2| , and the observed oscillations 



might be identifi ed as the Friedel oscillations. Also, with 
the above assumption, the density of states for electron 
tunneling into the compressible edge would be similar to 
the one for tunneling into a Fermi liquid (metal) . This is 
consistent with our intuitive expectations given the com- 
pressible nature of the system, if, loosely speaking, the 
characteristic energies for the motion of the charge and 
neutral (Fermi) part are comparable. 

We believe that it would be possible to construct an 
effective 1 + 1 dimensional theory along the edge, which 
has the same correlations that we expect, taking the co- 
ordinate normal to the edge to corresponds to time i.e. 
x = vt, and translating our diagramatic calculations into 
an effective interaction between a neutral and charged 
part. This would yield a model for the suppression of 
the correlation of the chiral boson theory (charge 
part), which assumes that its neutral and charge com- 
ponents move with the same velocity (v c = v n = v) 
along the edge. If the model is generalized to the one 
for which v c 3> v n , at sufficiently large momenta-high 
energies where the exchange part of the interaction is 
suppressed (due to a reduced overlap of the two one- 
dimensional spheres), we expect that the chiral boson 
correlations will be released. Therefore, the difference in 
the dynamics of the charge and neutral part appear to be 
a necessary condition for the decoupling of the edge and 
bulk (the charge and neutral part) at high enough ener- 
gies, as seen in experiments [|10f . (For a similar explana- 
tion of the simultaneous suppression of the neutral part 
see |12[| .) The "true" (low-energy) correlations should 
reflect the compressible nature of the system. 

In contrast with the edge problem, the bulk correla- 
tions in the compressible case that we considered seem 
to be similar to the ones in the incompressible case. The 
GM correlations are almost identical, and, due to the fi- 
nite screening, the quasiholes (correlation holes) have a 
chance to be considered as well-defined (albeit very ex- 
tended) objects (like skyrmions when the compressible 
degree of freedom - spin - is included in the incompress- 
ible problem p^|). This, intuitively, gives an additional 
support to the quasiparticle pictures of the bulk that 
we have by now. On the other hand, the edge correla- 
tions differ completely from the ones in the incompress- 
ible case. In incompressible states the edge physics is a 
reflection of the bulk physics, and the same quasiparti- 
cle picture of the bulk is possible on the edge. In the 
compressible case, and, in the plasma analogydue to the 
very weak screening on the edge, we probably can not talk 
about existence of the correlation hole which, in Read's 
picture of the bulk, attracts an electron, and creates a 
weakly-interacting composite object - a Fermi quasipar- 
ticle. In the scope of our approach, and in the first ap- 
proximation, electrons are unbounded and the edge of 
the compressible state appears to be similar to the edge 
of free electron gas (with an appropriate boundary con- 
dition). 



10 



ACKNOWLEDGMENTS 

We gratefully acknowledge useful discussions with A. 
Auerbach, J. Feinberg, J. H. Han, A. MacDonald, R. Ra- 
jaraman, S. Sondhi, A. Stern and X.-G. Wen, and espe- 
cially with N. Read. This work was partly supported by 
grant no. 96-00294 from the United States-Israel Bina- 
tional Science Foundation (BSF), Jerusalem, Israel, and 
the Technion - Haifa University Collaborative Research 
Foundation. M. M. also acknowledges support from the 
Fund for Promotion of Research at Technion, and the 
Israeli Academy of Sciences. E. S. would like to thanks 
the hospitality and support of the Institute of Theoreti- 
cal Physics in UCSB, Santa Barbara, where part of this 
work has been carried out. 



APPENDIX A: 

We consider the integral equation 
a general vertex 



for the case of 



•{- 



a-K5(qi - q f )} 



Then, the integral equation can be rewritten as 

7ra c b . 



(Al) 



V(qi,q f ) [1 
1 



b / dk 



i{qi - qf 

V(k,q f ) 



na 6(qi 



qf) 



i(k - qi )(k 2 + q 2 y 



(A2) 



where b = —Aum. If we try simply to iterate the equation 
in the limit when qi — > qf we find that each iteration 
produces a solution of the form 

V(qi,q f ) = 
at(qi,qf,q) 



-i(9i - qf) 

- q f ) + f(qi,qf,q) 



(A3) 



where a and p are fixed, i.e. do not change after some 
itrations, and f(qi,qf',q), keeps changing, but does not 
have any (new) singularity as qi — > qf. It can be assumed, 
from the iteration analysis, that f(qi, qf, q) is analytic in 
all of its variables in the long distance limit. 



If we assume that the solution is of the form A3 , and 



that a is an analytic function of its variables, the inte- 
gration on the right hand side can be done, and yields a 
new expression 



c{- 



1 



1 



(3b 



b n 



-i(qi-q f ) q 2 + q 

■f(Qi,qf) , f(iq,qf) i 



+ 7ra 8(qi - 



iq (iq 



b TT 
b 7T 



a (ft) 



«((?/) 



-i(qi - qf) qf + q 2 qj + q 2 

a(iq) 



q(iq-qi)(iq- qf) 



} 



(A4) 



where a(k) = a(k, qf; q). In order to equate the 6- 
functions on both sides, (at qi = qf), (3 must be 







(A5) 



q 2 +q'i - 



Then we equate the coefficients with — ^ r at the 

point qi = qf, to get 



*(qi,qf,q)\qi=q f 



c{q 2 + qjf 



{q 2 + qf — 7r a b c) 2 



(A6) 



To be consistent with the iteration result (and also with 
symmetry arguments), 



a(qi,qf,q) 



2 i 2 

q + qf 



q 2 + qj 



(q 2 + qf — 7r a b c) (q 2 + q 2 — tt a b c) 



(A7) 



although it is not consistent with our assumption that a 
is an an alyt ic function at the begining of the substitution 
of Eq. (A3). Still, it does satisfy the assumption in its 
long distance version 



a(g<, qf, q) ~ (?r - b c)2 (q 2 + q})(q 2 + qf) 



(A8) 



and the same approximation must be employed in the 
previous equations. 

To complete the solution we must find f{qi,qf',q) from 
the remaining equation 

[1 - n(a + 1) b c _L- 5 ]/(g i , q f ) = tt c b[ qf) i 



iq(iq - q z ) 

(A9) 



(We used a(iq) = which is consistent with the fact 
that the poles at k — ±iq in fc °^a a t the begining of 
the calculation were spurious.) In the long-distance (or 
small-momentum) approximation, i.e. when 



(o + l) 



f(qi,qf,q) _ f(iq,qf,q) 

(iq + qi) iq 



(A10) 



a nontrivial (nonzero) solution exists only when a = 1 
(Laughlin case). It is 



f(qi,qf,q) =q% + n 



(All) 



in the limit when qi qf (irrespective from the value of 

By power counting or by explicit calculation we can 
find out that the leading contribution in the Laughlin 
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case comes from the latter part of the solution. An in- 
troduction of an infrared cut-off is necessary, but the de- 
pendence on it disappears when the 8 - function part of 
the solution is included. When substituted in Eq. (35), 
this yields 



mln(£A) 



(A12) 



as the electrostatic energy of the impurity and its image 
counterpart in the longdistance approximation, where A 
is an ultraviolet cut-off (corresponding, e.g., to the in- 
verse screening length of the plasma) . There is no depen- 
dence on the infrared cut-off, because we are considering 
a half-plane (semi-infinite) system, and are recovering 
the well known result for that case. 



APPENDIX B: "IMAGE CHARGE" 
INTERACTION ENERGY IN A MODIFIED 
PLASMA 



where TZ denotes the real part, and r 2 = x 2 + y 2 . Note 
that this charge distribution decays algebraically towards 
the interior of the plasma, indicating its anomalously 
poor screening properties. The electrostatic energy as- 
sociated with the interaction of the impurity and the in- 
duced charge is then found to be (to leading order in 
small £) 



E„ 



m 2 A . A 
In — 

4<?o q c 



(B7) 



The higher order corrections decrease as a function of £. 
Multiplying by the inverse temperature = 2/m, and 
with the appropriate definition of the ultraviolet cutoff 
A, this result coincides with Eq. (Q), and hence is con- 
sistent with our diagramatic approach. 



We consider a point-like impurity of charge m, placed 
at a distance £ from the edge of a two-dimensional mod- 
ified plasma which occupies the half-plane x < 0. The 
plasma is characterized by a wave-vector dependent di- 
electric constant of the form 



<q) = i 



qo 
<1 



where qq = 



3 k 

4 7T 2 



I 



(Bl) 



and q = \q\ (see Eqs. ( [19 ) and (53)). The electrostatic 
potential generated by the charge distribution, that the 
external impurity induces in the plasma, is given by 

V lnd (r) = V sc (f) - V ex {r) ; (B2) 

here V sc (r) is the screened potential of the impurity 

V sc (r) — m / d 2 r'D(r, r') In \r — £x\ , where 



D(r,p) = 9{-x)0(-x') 



d 2 g 



-iq-{r — r') 



e-\q), (B3) 



and e(q) is given by Eq. (Bl) (the Theta functions re- 
strict the screening to the half-plane occupied by the 
plasma). This yields (in g-space) 



Vind{q) = - 



mnqoe 



Mi 



(q + qo)\q y \(\q y \ -n*) 



(B4) 



The (two-dimensional) Poisson equation then relates this 
component of the potential to the induced charge 



pind(q) = 7T V incl{q) 
Zir 



mq [\q y \ + iq x )e 



-\lyM 



(q + qo)\q y \ 

I) yields 



The Fourier transform of Eq. 



Ptod(f) = -0(-x)- TZ(r 2 -e + 2ifr) 



-3/2 



(B5) 



(B6) 
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